Pseudopotentials
Pseudopotentials in the Kleynman Bylander form¶
This page reports basic definitions and equations needed to evalute the nonlocal part of the Hamiltonian. We mainly focus on norm-conserving pseudopotentials. Providing a consistent introduction to the theory of pseudopotentials is beyond the purposes of this section, a more complete discussion of the topic can be found in specialized articles available in the literature [Hamann1979], [Bachelet1982], [Kleinman1982], [Hamann1989], [Troullier1991] [Gonze1991], [Fuchs1999]. For the generalization of norm conservation to multiple projectors (Optimized norm-conserving Vanderbilt pseudopotentials) we refer the reader to [Hamann2013]. The generation and validation of ONCVPSP pseudopotentials is described in the PseudoDojo paper [Setten2018].
For our purposes, we can summarize by saying that a pseudopotential in constructed in order to replace the atomic all-electron potential such that core states are eliminated and valence electrons are described by pseudo wavefunctions whose representation in the Fourier space decays rapidly.
Modern norm-conserving pseudopotentials are constructed following the procedure described in [Hamann1979], [Hamann1989], [Troullier1991] such that the scattering properties of the all-electron atom are reproduced around the reference energy configuration, up to first order in energy. In modern ab initio codes, the fully separable form proposed by Kleynman and Bylander is usually employed as it drastically reduces the number of operations required for the application of the nonlocal part of the Hamiltonian as well as the memory required to store the operator [Kleinman1982]. In the KB form, the interaction between a valence electron and an ion of type is described by means of the operator:
where is a purely local potential with a Coulomb tail .
Note
For simplicity, the discussion is limited to pseudopotentials with a single projector per angular channel. The generalization to multi-projector pseudopotentials requires introducing an additional index in the KB energies i.e. and in the projectors .
The so-called Kleinman-Bylander energies, , measure the strength of the nonlocal component with respect to the local part. The projectors are short-ranged functions, expressed in terms of a complex spherical harmonic multiplied by an - and atom-dependent radial function
The total nonlocal part of the Hamiltonian is obtained by summing the different atom-centered contributions. The final expression in real space reads:
where runs over the sites of the Bravais lattice, and over the positions of the atoms with the same type located inside the unit cell. Due to the nonlocality of the operator, its Fourier transform depends on two separate indices, and , instead of the simple difference . The shorthand notation is used in the following to simplify the derivation.
Using the Rayleigh expansion of planewaves in terms of spherical harmonics and spherical Bessel functions :
the Fourier representation of the projector defined by \ref{eq:KB_projectors} can be expressed as:
Finally, the expression for the total nonlocal operator in reciprocal space is:
Note
Abinit employs iterative eigenvalue solvers to solve the KS equations. This means that we only need to compute i.e. we only need to apply the Hamiltonian onto a set of trial eigenvectors. Therefore the full matrix is never constructed explicitly.